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Abstract 

Recently it has been shown by the present author [H. Hasegawa, Phys. 
(— I , Rev. E (in press): larXiv:0904.2 399j that the interpolation approximation 

O \ (lA) to the generalized Bose-Einstein and Femi-Dirac distributions yields re- 

sults in agreement with the exact ones within the 0{q — 1) and in high- and 
low-temperature limits, where ((/—I) expresses the non-extensivity: the case of 
g = 1 corresponding to the conventional quantal distributions. In this study, 
we have applied the generalized distributions in the lA to typical nonextensive 
subjects: (1) the black-body radiation, (2) the Bose-Einstein condensation, 
(3) the BCS superconductivity and (4) itinerant-electron (metallic) ferromag- 
netism. Effects of the non-extensivity on physical quantities in these nonex- 

■ tenisive quantum systems have been investigated. A critical comparison is 
^ , made between results calculated by the lA and the factorization approxima- 

I tion (FA) which has been so far applied to many nonextensive systems. It 

has been pointed out that the FA overestimates the non-extensivity and that 
, it leads to an inappropriate results for fermion systems like the subjects (3) 

> • and (4). 
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1 Introduction 



Considerable works have been made on the nonextensive statistics since Tsalhs 
proposed the generahzed entropy (called the Tsallis entropy) |1] which is a one- 
parameter generalization of the Boltzmann-Gibbs entropy with the entropic index 
q: the Tsallis entropy in the limit of g = 1.0 reduces to the Boltzmann-Gibbs entropy 
(for a recent review, see [2J). In recent years, much attention has been paid to an ap- 
plication of the nonextensive statistics to quantum subjects, in which the generalized 
Bose-Einstein and Fermi-Dirac distributions (called g-BED and g-FDD hereafter) 
play important roles. The four methods have been proposed for g-BED and g-FDD: 
(i) the asymptotic approach (AA) [3j obtained the canonical partition function valid 
within 0((g — 1)//cbT), (ii) the factorization approach (FA) [Ij employed the de- 
coupling, factorization approximation in evaluating the grand-canonical partition 
function, (iii) the exact approach (EA) [5l E] derived the formally exact expression 
for the grand canonical partition function expressed in terms of the Boltzmann- 
Gibbs counterpart, and (iv) the interpolation approximation (lA) [7] was proposed 
based on the EA, yielding results in agreement with those obtained by the EA within 
0{q — 1) and in high- and low-temperature limits. Among the four methods, the 
FA has been mostly adopted in many quantum subjects including the black-body 
radiation [SI El [ID!, early universe [HI [12], the Bose-Einstein condensation (BEG) 
[13]- [20], metals [21], superconductivity [221 [23], spin systems [24J-[29J and itinerant- 
electron (metallic) ferromagnets [30j . This is due to a simplicity of the expression 
of the generalized quantal distributions in the FA. 

Quite recently, however, it has been pointed out that the FA is not accurate from 
a study of the EA [7] . It might be necessary to examine calculations previously made 
with the use of the FA by a new calculation with the lA, which is the purpose of 
the present paper. We will discuss the four quantum subjects: (1) the black-body 
radiation, (2) the Bose-Einstein condensation, (3) the BGS superconductivity and 
(4) the itinerant-electron ferromagnets, to which the FA has been applied [H [9l 
[in][I3]-[2U] [221E3][3n]. 

The paper is organized as follows. In Sec. 2, we briefly discuss the g-BED and 
g-FDD in the EA and lA after Ref. [7|. Then we apply the lA to the four sub- 
jects mentioned above. The black-body radiation is discussed in Sec. 3, where the 
g-dependent Stefan-Boltzmann coefficient and the Wien shift law are calculated. In 
Sec. 4, we investigate the BEG, calculating the critical temperature and the temper- 
ature dependence of the energy and specific heat as a function of g. The g-FDD in 
the lA is applied to the BGS superconductivity in Sec. 5, where the g dependence 
of the superconducting critical temperatures and the characteristic ratios related 
with the ground-state order parameter, critical temperature and specific heat in the 
BGS theory. In Sec. 6, we discuss magnetic and thermodynamical properties of the 
itinerant-electron ferromagnets described by the Hubbard model combined with the 
Hartree-Fock approximation. In Sec. 7, we present qualitative discussions with the 
use of the generalized Sommerfeld low-temperature expansion. Section 7 is devoted 
to our conclusion. 
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2 Generalized quantal distributions 
2.1 An exact approach 

When the maximum entropy method with the optimal Lagrange multipher [31] is 
adopted, the generahzed quantal distribution for the state k (whose number operator 
is hk) is given by [7] 

fg{ek,P) = -^Tr{[l + {q-l)P{H-fiN-E^ + fiN^)]T^^hk}, (1) 

for g > 1, (2) 
for g < 1, (3) 



27rX 



with 



A', = rr{|l + (,-l)/3(ff-^iV-£;, + ,,JV,)]T^}, (4) 
= rG(«;^,|^)e~)E.M.„ fo.,>l, (5) 



where 



k 

ni{u) = i-V ln[lTe-"(^^-'^)], (8) 

k 

G(«;a,6) = -^w--ie-^", (10) 
r(a) 

H{t;a,b) = T{a + l)b-'' i-ty-^e-^'. (11) 

Upper and lower signs in Eqs. ([7])- ([9]) are applied to boson and fermion, respectively, 
r(z) stands for the gamma function and C denotes the Hankel path in the complex 
plane [3, E]- It is noted that Eqs. ([I])-([SD include Ng and Eg which should be 
determined in a self-consistent way [see Eqs. (2) and (3) in Ref. 0]. Such self- 
consistent calculations have been reported for the band electron model and the 
Debye phonon model in Ref. 
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2.2 The interpolation approximation 

Self-consistent calculations for fq{ek, P) including Ng and Eg are rather tedious and 
difficult . In order to overcome this difficulty, we have proposed the lA |[7j , assuming 
that 

^e"(^^-^^«)SiH = 1, (12) 



in Eqs. ([2]) and (jSj). Then the generalized distribution in the lA is given by 
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[{l-q)(3]-—^ V27ry Jc 



for q > 1.0, (13) 

(_t)-T^e~'ra3 /i(e, -t) dt 

for g < 1.0. (14) 



g-BED 

With the use of Eqs. f|T3|) and f|T4l) . the analytic expression of the g-BED in the 
lA is given by [7j 



oo 



/f(e,/?) = J^K'^^+^^I" for0<g<3, (15) 



n=0 

where expresses the g-exponential function defined by 

Cg = expg(x) = [1 + (1 - g)x] 1^ for 1 + (1 - g)x > 0, (16) 
= for 1 + (1 - q)x < 0, (17) 

with the cut-off properties. 
g-FDD 

Similarly, the analytic expression of the g-FDD in the lA is given by [7j 

//^(e,/?) = F(e,/3) for e > /i, (18) 

= I fore = /x, (19) 

= 1.0- F{\e- fi\+ fi,(3) ioreKfi, (20) 

with 

oo 

F(e,/3) = ^ (-l)«[e-("+i)-]'' for < g <3. (21) 

n=Q 

Note that = [1 - (1 - g)(n + l)x] ^ [e-^]("+i) in Eqs. and ([2lD 

of g ^ 1.0 where — >■ 



except for g = 1.0. fg {e, (3) given by Eqs. (fT5l) - (!2TI) reduces to fi{e, (3) in the limit 



Table 1: Generalized quantal distributions in the limits of g — >^ 1, /3 — ^ cxd and /3 — 



method 




(3^oo (FDD) 


/5 -> 


EA" 


fi + {q- 1) 






e(/i-6) 




lA^ 


fi + {q- 1) 






e(/i-6) 




FA'^ 


A-i(g-l)/3(6-/i)2^ 


e(/i-e) 


fiq 


FAq'^ 


/i + (g-l)[(6-/i)-|/3(e-/in^ 


e(/i-6) 





/i = l/(e'^'''^~'^'' =F 1): 0(x), the Heaviside function: e^, g-exponential function. 
" the exact approach [7] 

the interpolation approximation [7] 

the factorization approximation [4] 

the factorization approximation [Eq. ( I23l) ] 



On the contrary, the g-BED and g-FDD in the FA are given by [1] 

^ oo 

fn^^^) = ^ 1 = (22) 

where the upper (lower) sign is applied to boson (fermion). It is noted that if we 
adopt a factorization approximation: e^''"^^''^ ~ [e"^]'-""'"^^ in Eqs. f ll5l) - fl2Tl) . we 
obtain 

oo ^ 

/f^^(e,/3) = T^T^pT^' (23) 

which is similar to Eq. (1221) and which is referred to as the FAq hereafter. 

A comparison among the 0{q — 1) contributions to the generalized quantal dis- 
tributions in the EA, lA, FA and FAq is made in Table 1. It is stressed that the 
lA has the interpolation character yielding good results in the limits of g — > 1.0, 
/9 — > cxD and (3 — > 0.0 [7]. The 0(g — 1) contributions in the FA and FAq are different 
from that in the EA. In the limit of /3 ^ 0.0, f^^ie), f^^{e) and /^^^^(e) reduce to 
[e"^^]"^, while f^^ reduces to e"'^'^. In the limit of /3 — > oo, all the g-FDD become 
0(/i — e), where 9(x) denotes the Heaviside function. More detailed comparisons 
among various methods have been made in Ref.[7]. 
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3 Black-body radiation 



We will apply the g-BED given by Eq. ( flSl) in the lA to the black-body radiation 
model with the photon density of states per volume given by 

p(cu) = Cu?, (24) 

where C = l/vr^c^ and c denotes the light velocity. The generalized Planck law is 
given by 

Dgicu) = hjp{u;)f^^{hj,P). (25) 

The g-BED in the FA was adopted to the black-body radiation in Refs. [HI M, [ID] • 

The g-BEDs (with /i = 0.0) calculated in the lA and FA are shown by solid 
and dashed curves, respectively, in Fig. [It^a) with the logarithmic ordinate: they 
are indistinguishable in the linear scale. For q = 1.2, tails of g-BED obey the 
power law. In contrast, for q = 0.8, g-BED has a compact form with the cut-off 
behavior: fq{hjj) = 0.0 for pkuj > 5.0. Solid and dashed curve in Fig. [D^b) express 
the generalized Planck law Dq{uj) calculated by the lA and FA, respectively. For 
q = 1.2, the distribution of Dqiuj) in the high-frequency region is much increased. 
This trend is reversed for q = 0.8. The effect of the non-extensivity in the FA is 
much overestimated compared to that in the lA. 
We obtain the generalized Stefan-Boltzmann law. 



Eq = / Dq{uj)duj, (26) 
Jo 

= (27) 



with 

r(-i, - 3) 



(g-i)^r(^ + i) 
r(^ + i) 



(i-^)'r(T^ + 5) 
1 

(2-g)(3-2g)(4-3g) 



for q > 1.0, (28) 

for q > 1.0, (29) 
for < g < 4/3, (30) 



where cxi is the Stefan-Boltzmann constant for g = 1.0. The g dependence of cTg 
calculated in the lA and FA are shown by solid and dashed curves, respectively, 
in Fig. [2] with the logarithmic ordinate. With increasing g, a is monotonously 
increased. 

Substituting Eq. f[T^ to Eq. fl25]) . we obtain Um where Dq{ijj,(3) has the maxi- 
mum. 
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Ol v^oo 1 r -(n+l)/3?k<Jmi„ 



^ y Z^n=0 n! l^'} J 

for 1.0, (33) 

whose solution expresses the generahzed Wien shift law. The solid curve in Fig. 
[3] shows the calculated ratio of ujm,q/^m,i as a function of q. With increasing q 
above q = 1.0, the ratio is increased whereas it is decreased with decreasing q below 
unity. Dashed and chain curves show the results of the FA and AA [uJrn.q/^m,i = 
1 + 6.16 (g — 1)] [3], respectively. 

4 Bose-Einstein condensation 
4.1 Basic equation 

We will study this subject, applying the lA to a bose gas with the density of states 
given by 

p(e) = A6^ (34) 

where r = d/2 — 1 for ci-dimensional ideal bose gase, r = d — 1 for bose gase trapped 
in ci-dimensional harmonic potential, and A stands for a relevant coefficient. The 
nonextensive BEC was discussed in Refs. fl3]-[20] with the use of the FA. 

By using the g-BED given by Eq. (fT5|) in the lA, we obtain the number of 
electrons given by 



with 



N = N, + N„ (35) 



— = J2 e~^"+^^", for g = 1, (36) 



e 

n=0 



oo 



= ^[e-("+^)"]' for0<g<3, (37) 

n=0 

Ar(r + i)r(^ -r) 

= {ksTr' _ ^y^,^^^ ^ Mr + 1, a) for 1< g < 3, (38) 
= {kBTy+^AT{r + l)(f){r + l,a) for g = 1, (39) 

^ Ar(r + i)r(^ + i) 

= iksTr' (1 _ ,).+ir(^ + r + 2) '^^^^ + 1- ^) < ^ < 1- (40) 



where a = —jSfi (> 0), and Nc and A^e denote the numbers of electrons in the 
condensed and excited states, respectively. Here 0q(z, a) is the generalized Bose 
integral defined by 



oo !^-na]z-{z-l)q 



{z,a) ^ JZ^^^i forgf?^>l, (41) 



n 

n=l 
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which reduces to 



oo ^ 

(j)q{z,a) 5Z;^ = C(^) for a ^0.0, 



n=l 

oo 

E 

n=l 



for q 1.0, 



(42) 



(43) 



^(2;) and (f){z, a) being the Riemann zeta function and the Bose integral, respectively. 
4.2 Critical temperature 

The number of electrons in the excited state is bounded by Eqs. f l38|) - fH0|) with 
a = 0.0. Then the critical temperature of the BEC, Tc, below which a vanishes is 
given by 



A^r(-i^ + i) 



Ar(r + l)C(r + l)r(^- 
N 



_v4r(r + l)C(r + 1)_ 



for q = 1, 



iVr(^ + r + 2) 



Ar(r + l)C(r + l)r(^ + l) 



for 1 < g < 3, (44) 
(45) 

for < g < 1, (46) 



We note that = for r = {i.e., free boson with d = 2 or d = 1 boson with 
harmonic-potential) because ({1) = 00. Equations ( l4^ - (H6l) lead to 



T 



r(^ + r + 2) 

r(T^ + i) 



r+l 



for 1 < g < 3, 



for < g < 1. 



(47) 



(4^ 



The solid curve in Figs. H] (a) and (b) show the g dependence of the ratio of Tc^q/Tc^i 
for r = 1/2 and r = 2, respectively, calculated by Eqs. ( H7I) and ( l48l) . The critical 
temperature is decreased with increasing g. 

On the other hand, the critical temperature in the FA is given by [T7] 



T 



gr + i) 

C,(r + 1) 



with 



Ur + 1] 



C(r + 1 
1 



1 ^ r(^-r-i) 

' n=l "-yq-l! 



r(^ + r + 2) 



for 1 < g < 3, 
for g = 1, 
for < g < 1. 



(49) 

(50) 
(51) 
(52) 
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Dashed curves in Figs. H] (a) and (b) express the results of the FA calculated by 
Eqs. ( H9|) - (!52|) . The effect of the non-extensivity is overestimated in the FA: T^^ 
vanishes at g > 1.67 and q > 1.33 for r = 0.5 and r = 2.0, respectively. In contrast, 
T/"^ vanishes at g > 3.0 and g > 1.5 for r = 0.5 and r = 2.0, respectively. 

4.3 Condensed states 

The number of electrons in the condensed states A^^^ is given by 



which depends on r but it is independent of q. 

4.4 Energy and specific heat 

The total energy is given by 

v4r(r + 2)r(^ - r - 1) 
E = iksTY^' (q - ly+^ fi^ + 1) '^^(^ + ^'") forg^l, (54) 

= (fcBT)"+2^r(r + 2)0(r + 2,e-") for g = 1. (55) 

Above Tc, a is temperature dependent because it is adjusted as to conserve the 
total number of electrons. 

Then its temperature dependence is given by 

da (r + 1) (j)q{r + l,a) 



dT (r + 1 — rg) T 0g(r, a) 
(r + 1) 0(r + l,a) 
T (hir.a) 



for g 7^ 1, (57) 
for g = 1. (58) 



Taking into account the temperature dependence of a, we obtain the specific heat 
at T > Tc given by 

C [TV^^ (r + 1) 

- ' ^ ^ ^ {(r + 2)0,(r + 2,a) 



keN \TJ [r + 2-(r + l)g]C(r+l) 

(r + l)[r + 2- (r + l)g][0g(r + l,a)]^ 
0g(r,a) 



} for g ^ 1, (59) 



T\'+' (r + 1) r , , , , J0(r + l,a)12 

^ ^ ' (r + 2)0(r + 2,a)-(r + l)L^^ ' 



Tc) C(r + 1) 



0(r,a) 
for g = 1. (60) 
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Below Tc where a = 0.0, we obtain the specific heat given by 
C (^ + i)(^ + 2)C(r + 2) 



keN \TJ [r + 2-(r + l)g]C(r + l^ 

TV^' (r + l)(r + 2)C(r + 2) 



Tc) C(r + 1) 



for q^l, (61) 
for q = l. (62) 



The calculated specific heats for r = 0.5 and 2.0 are plotted in Figs. O^a) and 
(b), respectively. The magnitude of the specific heat is monotonously increased with 
increasing q and/or r. A jump in the specific heat at is given by 

AC _ C{T, - 0) - C{T, + 0) 

(r + lYCir + 1) 
= ^ I ' for0<g<3. (64) 

C(r) 

Equation fl64l) shows that for r < 1.0, AC vanishes and C is continuous at because 
of the divergence in C{r). Then AC vanishes for r = 0.5 in Fig. [T2l(a) while it is 
finite for r = 2.0 in Fig. lT2\h). 



5 BCS supercondunctivity 

5.1 Model hamiltonian 

The BCS model hamiltonian is given by 

H = ^ek{nk + n^k) + '^Vkk'bl,bk, (65) 

k kk' 

where the operator bl (= c|.cLfc) creates a Cooper pair in a singlet superconducting 
state. The attractive interaction Vkk', which is the origin of the superconductivity, 
is assumed to be effective only near the fermi level fi, 

Vkk' = -V for |efc - |efc/ - yu| < (66) 
= otherwise, (67) 

where Iiwd = ksQu, and ujd and Qd denote the Debye energy and temperature, 
respectively. The BCS model was discussed with the use of the FA [22l[23]. We will 
examine nonextensive superconductivity, by using the g-FDD in the lA. 



5.2 Order parameter 

By using the double-time Green function method for the nonextensive quantum 
systems [5l [6], Numes and Mello [22] derived the self-consistent equation for the 
order parameter A given by 

A = Vp [f^{-E,P)-f,{E,(3)]de, (68) 
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with 



E = + A2, 



(69) 



where p = p(/i) expresses the density of states at the fermi level. The ground-state 
order parameter Aq is determined by 



1 



de. 



(70) 



Because the g-FDD at T = is the same as /i(e) 0, the ground-state order param- 
eter is identical with the BCS result for g = 1.0: 



Ao 



-l/Vp 



(71) 



independently of q. The critical temperature 
given by 



1/kBPc) where A vanishes is 



1 



2/5(6, /5c 



de. 



(72) 



We show in Fig. [6l the g-FDD calculated by the lA and FA for q = 0.8 and 1.2, 
the result for q = 1.0 being plotted by the chain curve. For q = 1.2, a tail of the 
distribution at large e obeys the power law. In contrast, the distribution for q = 0.8 
has a compact form with the cut-off properties: /q(e) = 1.0 for /?(e — p) < —5.0 
ad fq[e) = 0.0 for [3{e — p) > 5.0. These properties in the g-FDD distribution are 
more clearly realized in its derivative, —dfq{e)/de, which is plotted in the inset of 
Fig. [61 We note that —dfg{e)/de in the lA is symmetric with respect to e = p 
independently of q, while that in the FA is not for q 7^ 1.0. 

We have solved the self-consistent equation fl68p . by using the g-FDD in the lA. 
Figure [7|(a) shows the temperature dependence of the order parameter for various q 
values with Vp = 0.3. With decreasing q from unity, the temperature dependence 
of the order parameter become significant and the Curie temperature is decreased. 
On the contrary, with increasing q above unity, we observe the opposite tendency: 
the temperature dependence of the order parameter become less significant and the 
Curie temperature is increased. The properties of the order parameter are clearly 
realized in Fig. Wljo), where the normalized order parameter is plotted. The q 
dependence of the normalized critical temperature, Tc^q/Tc^i, is plotted in Fig. [HI 
showing an monotonous increase in with increasing q. The dashed curve in Fig. 
[H] shows the q dependence of the critical temperature calculated with the g-FDD in 
the FA, which shows the different q dependence of Tc from that in the lA. The chain 
curve will be discussed shortly. 

We may obtain the (q — l) expansion of the superconducting critical temperature, 
by using the {q — 1) expansion of /f^(e) given in Table 1, leading to 



2/, 



I A/ 



2/1 



X 



X 



2(g-l) 



dx 2 dx"^ 



+ 



(73) 
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Substituting Eq. (!73|) to Eq. (!72l) . we obtain 



from which the critical temperature is given by 

(75) 

Tc,i [={2e^ /tt) hojD e~^^^^ ] denoting the critical temperature for q = 1.0 and 7 
(=0.577) the Euler constant. The chain curve in Fig. [8] expresses Eq. f l75|) . which 
is in good agreement with the solid curve expressing the result obtained by a self- 
consistent calculation of Eq. (I68l) . 



c,q 



TcA 



1) + 



5.3 Specific heat 

The specific heat in the superconducting states is given by 

C, . 1 r - Me, (76) 



T Jo V 2 dT J \ dE 
while that in the normal states is given by 



Figure [HI shows the temperature dependence of the specific heat for various q val- 
ues calculated with the use of g-FDD in the lA. The temperature dependence of 
the specific heat shows the well-known exponential temperature dependence in the 
superconducting states for q = 1.0. When q is decreased from unity, this behavior 
becomes more significant. On the contrary, when q is increased from unity, the 
exponential temperature dependence is changed to nearly linear and then convex 
behavior. There is a jump in the specific heat at given by 



AC = C,(T,-0)-C,(T, + 0), (7 



dT Jq \ de 

It is interesting to investigate the q dependence of the ratios given by 

2Ao /AC\ 



(79) 



^0) 



which are 3.53 and 1.43, respectively, for q = 1.0 in the BCS theory. Solid curves 
in Fig. [To] express the q dependence of these ratios calculated in the lA, showing 
that both the ratios are deceased with increasing q. In contrast, the q dependence 
of the ratios calculated in the FA shown by dashed curves are rather different: with 
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Table 2: Ratios of 2Ao//c_bTc and t^CjCn of typical BES superconductors" and 
estimated q values with the use of Eq. flSTj) 



materials 


2Ao/A;BTe 




q 


(BCS) 


3.53 


1.43 


1.00 


Cd 


3.44 


1.32, 1.40 


1.05 


Zn 


3.44 


1.3 


1.05 


V 


3.50 


1.49 


1.01 


Al 


3.53 


1.29-1.59 


1.00 


Sn 


3.61 


1.60 


0.96 


Tl 


3.63 


1.5 


0.94 


Ta 


3.63 


1.59 


0.94 


In 


3.65 


1.73 


0.93 


Nb 


3.65 


1.87 


0.93 


Hg 


3.95 


2.37 


0.76 


Pb 


3.95 


2.71 


0.76 



Ref. [32] 



increasing g, AC/C„ is decreased but ^t^^jkBTc is increased. The q dependence in 
2^Q/kBTc is due to the g-dependent Tc shown in Fig. [8] because Aq is independent 
of q. 

Table 2 shows observed ratios of 2/S.q /kBTc and /S.C jCn of typical superconduct- 
ing materials [32]: 2Ao/A;bTc vs. /S.C jCn plot of these data is presented in Fig. [TTl 
We note the tendency such that materials with larger 2Ao/A;bTc have larger AC/C„, 
which is consistent with our calculation shown by the solid curve, except for Hg and 
Pb which are considered not to be weak-coupling superconductors. From Eq. (!75ll . 
we obtain 

Applying Eq. (IHT!) to the observed values of 2/S.q /kBTc, we have estimated q values of 
the materials, which are shown in the fourth column in Table 2. A similar estimate 
of q values for some BCS materials was made in [22j but using the FA. The FA, 
however, yields the tendency such that materials with larger 2Ao/fcBTc have smaller 
AC/Cni ill contrast with our result in the lA. 



(81) 
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6 Itinerant-electron ferromagnets 



6.1 The Hartree-Fock approximation 

We will discussed itinerant-electron (metallic) ferromagnets described by the Hub- 
bard model given by [30] [33] 

a i (7 i,j i i 

(82) 

Here rijo- = al^ai^, (af^.) denotes an annihilation (creation) operator of a a-spin 
electron (a =t, |) at the lattice site i, eo the intrinsic energy of atom, tij the electron 
hopping, U the intra-atomic electron-electron interaction, B an applied magnetic 
field and fiB the Bohr magneton. With the use of the Hartree-Fock approximation, 
Eq. (1821) becomes the effective one-electron Hamiltonian given by 

H = 5^ 5^ eo + J] 5^ tij alaj^ + U '^(Mrii^ + {ni^)nii) 

a i cr i,j i 

-fiB B'^im^ - mi), (83) 

i 

where the bracket (■) denotes the expectation value [Eq. fl86|) ]. This subject was 
previously discussed in Ref. [30] with the use of the FA. 



6.2 Magnetic moment 

Self-consistent equations for the magnetic moment (m) and the number of electrons 
{n) per lattice site are given by 



m = (ni)-(ni), (84) 

n = {ni) + {ni), (85) 

with 

M = I p.{e)fg{e)de, (86) 

PTa(e) = Po (^e - eo - -^(riTm) i/ifiS^ , (87) 

P^('^ = (88) 
k 

where /g(e) expresses the g-FDD given by Eqs. (fT8|) - (l2T|) . po(e) denotes the density 
of states, efc is the Fourier transform of tij and Na the number of lattice sites: the 
plus and minus signs in Eq. (187|) are applied to j- and |-spin electrons, respectively. 
From Eqs. (j84l) - (!88l) . m and fi are self-consistently determined as a function of T 
for given parameters of q, n and U and density of state, po(e)- 
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We have performed model calculations bearing in mind Fe, which has seven d 
electrons and the ground-state magnetic moment of 2.2 hb- By using a bell-shape 
density of states for a single band given by [30] 




we have adopted U/W = 1.75 and n = 1.4 electrons as in [30], where W denotes a 
half of the total bandwidth. 

Fig. [12] shows the temperature dependence of the magnetic moment m for 
q = 0.8, 1.0 and 1.2 calculated by the lA and FA. For q = 1.2, the temperature 
dependence of magnetic moments becomes more significant and the Curie temper- 
ature becomes lower than for q = 1.0 in the lA. On the other hand, for q = 0.8, the 
temperature dependence of m becomes less significant and the Curie temperature 
becomes higher than for q = 1.0 in the lA. The behavior of m in the FA is quite 
different from that in the lA: the Curie temperature is more decreased both for 
g = 0.8 and 1.2 than for q = 1.0. This fact is more clearly seen in Fig. [T31 where 
Tc is plotted as a function of q. The Curie temperature in the lA monotonously 
decreased with increasing q. On the contrary, Tc in the FA is almost symmetric 
with respect to g = 1.0 where we obtain the maximum value of ksTc/W = 0.143. 
If we adopt W ~ 2.5 eV obtained by the band-structure calculation for Fe [31], the 
calculated Curie temperature at g = 1.0 is Tc ~ 3500 K, while the observed Tc of 
Fe is 1044 K [351. 



6.3 Energy and Specific heat 

The energy per lattice site is expressed by 

E = Je[p,ie)+pMmde-^in'-m'), (90) 

from which the electronic specific heat is given by 

dE dE dE dm dE dfi , 

dT ^ df ^ d^'dT ^ ~dJIdT' ^ ^ 

with 

- / e(e-M)h(e)+Pi(6)]^d6, (92) 

[M^)-Pli^)]^de, (93) 

Ihi^) +PM^de. (94) 

Analytic expressions for dm/dT and dfi/dT in Eq. (!9T|) are given by Eqs. (A.3)- 
(A.8) in Ref. [30] . 

Fig. [T^ shows the temperature dependence of the specific heat C for g = 0.8, 1.0 
and 1.2, calculated by using the lA and FA. In the lA, C for g = 0.8 is smaller than 
that for g = 1.0. In contrast, C in the FA of g = 0.8 is larger than that of g = 1.0. 



dE 


1 




~T 


dE 


U 


dm 


~ 2" 


dE 




dp 


-/ 
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6.4 Spin susceptibility 

The spin susceptibility is expressed by [30 



dm , , 



from which the paramagnetic spin susceptibihty is given by 



2 l(\r\ 

^ = 71 JT^^ 96 

(1 - ^Xo) 



with 



Xo = -/ P(e)^rfe. (97) 

Figs. [15] shows the temperature dependence of the inversed susceptibihty 1/x 
calculated by the I A and FA for q = 0.8, 1.0 and 1.2. The Curie temperature Tc, 
which is realized at 1/x = 0, is monotonously decreased with increasing q in the lA, 
which is different from its q dependence in the FA, as shown in Fig. 



7 Discussion 

We will qualitatively elucidate the difference between the results calculated with 
the lA and FA for itinerant-electron ferromagnets. The generalized Sommerfeld 
expansion including an arbitrary function 0(e) and the g-FDD fq{e) is given by 



b{e)fg{e) de, (98) 
0(e) de + J2 Cn,, (A;BT)>('^~i)(/i), (99) 



n=l 



with 



.^/(e-,r^*, (100, 



which is valid at low temperatures. Expansion coefficients for q = 1.0 are given 
by C2,i = 7rV6 (=1.645), C4,i = 77rV360 (=1.894), and c„,i = 0.0 for odd n. The 
coefficients Cn,q for n = 2 and 4 in the lA are given by [7] 

= J^y (101) 

C2,i (2 - q) 

(2-g)(3-2g)(4-3g)' ^ ' 
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whereas c{"^ = = 0. Results in the lA is in agreement with those of the EA 
within 0{q-l) [t]. 

On the other hand, the FA yields [30| 



6 

7tt 



4 



60 



where a ~ 9.2, b ~ 110 and the 0{q — 1) contributions to cf^^ and cf^ are vanishing 
[30]. cf^ and cf^ are not zero, which is in contrast with the results of the EA and 
lA. This is due to a lack of the symmetry in —df^^{e)/de, as shown in the inset of 
Fig. El 

Figure [11] shows the q dependence of c„_q for n = 1 — 4 calculated by the lA and 
FA. We note that the q dependence of cf ^ and cf^ is symmetric with respect to 
g = 1.0 whereas that in the lA is not. 

By simple calculations using Eqs. (jH!]), ([85]), ([90]), ([96]), ([99]) and f[TOO]) . we obtain 
the magnetic moment m{T), the specific heat C at low temperatures and the Curie 
temperature Tc,g given by [30] 

m(T) ~ m{0)-aT^, (103) 

C{T) ^ 7,T, (104) 

with 

a = C2,q{p[-p\), (106) 

7, = 2c2,j2(p^+pJ-f/m(0)(pt-Pi)], (107) 

where = Pa{p), P = pip), P^^"* = (Pp{p)/de'^, and m(0) is the ground-state 
magnetic moment. Equation (11051) leads to 



rr / \ -1/2 



TcA \C2A 



(108) 



~ 1 - -^^-^-^ + ■ ■ in the lA, (109) 
nin - 1^2 

~ + in the FA. (110) 

Equation (1109^ shows that with increasing (decreasing) q from unity, the Curie 
temperature in the lA is decreased (increased). In contrast, Eq. (11 101) shows that 
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the Curie temperature in the FA is decreased with increasing \q — 1|. These are 
consistent with the results shown in Fig. [131 The coefficient of C2,g expresses the 
contribution from the Stoner excitations, which play important roles in magnetic and 
thermodynamical properties of itinerant-electron ferromagnets. It is noted from Eqs. 
fl3UI) and fll02p that the Stefan-Boltzmann constant cXg is related to the Sommerfeld 
expansion coefficient c^^g as given by ag/ai = 04^5/04^1. Thus the difference in the 
Sommerfeld expansion coefficients in the lA and FA reflects on the difference in the q 
dependence of the physical quantities calculated by the two kinds of approximations. 

8 Conclusion 

By using the generalized distributions of g-BED and g-FDD in the lA which was 
proposed in Ref. we have discussed four nonextensive quantum subjects: (1) the 
black-body radiation, (2) Bose-Einstein condensation, (3) the BCS superconductiv- 
ity and (4) itinerant-electron ferromagnetism. The effect of the nonextensivity of 
(g — 1) has been shown to be appreciable in these systems. A comparison between 
the results obtained by the lA and FA has shown that the FA generally overesti- 
mates the effect of the nonextesivity of Ig — 1| and that the g-FDD in the FA may 
yield qualitatively inaccurate results for fermion systems such as (3) and (4). This 
is considered to be due to the inappropriate g-FDD in the FA. It would be necessary 
to examine the validity of the previous studies which have made with the use of the 
FA. 
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Figure 1: (Color online) The uj dependence of (a) fq{huj) and (b) Dq{uj) calculated 
by the lA (solid curves) and FA (dashed curves) for q = 0.8 and 1.2, result for 
g = 1.0 being plotted by chain curves. 



Figure 2: (Color online) The q dependence of the coefficient of the generalized 
Stefan-Boltzmann law, cr^, calculated in the lA (the solid curve) and FA (the dashed 
curve). 



Figure 3: (Color online) The q dependence of the generalized Wien shift law, Um,q, 
calculated in the lA (the solid curve), FA (the dashed curve) and AA (the chain 
curve) ^. 



Figure 4: (Color online) The q dependence of the critical temperature of the Bose- 
Einstein condensation, Tc^g, for (a) r = 0.5 and (b) r = 2.0 calculated by the lA 
(solid curves) and FA (dashed curves). 



Figure 5: (Color online) The temperature dependence of the specific heat C for (a) 
r = 0.5 and (b) r = 2.0 calculated in the lA, Tc denoting the critical temperature 
for the EEC. 



Figure 6: (Color online) The e dependences of the g-FDDs of fq{e) [—dfq{e)/de 
in the inset] for q = 0.8 (solid curves) and 1.2 (bold solid curves) in the lA, and 
those for q = 0.8 (dashed curves) and 1.2 (bold dashed curves) in the FA, results 
for g = 1.0 being plotted by chain curves. 



Figure 7: (Color online) (a) The temperature dependence of the superconducting 
order parameter A, and (b) the temperature dependence of the normalized order 
parameter A calculated by the lA for various q values (Vp = 0.3). 



Figure 8: (Color online) The q dependence of the critical temperature Tc^q/Tc^i 
calculated by the lA (the solid curve) and FA (the dashed curve) with V p = 0.3, 
the chain curve denoting the result of (g — 1) expansion given by Eq. ( l75l) . 
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Figure 9: (Color online) The temperature dependence of the specific heat A calcu- 
lated by the lA for various q values {V p — 0.3). 



Figure 10: (Color online) The q dependence of the ratios of 2Ao//cbT'c and 
AC/Cn{Tc) calculated by the lA (solid curves) and FA (dashed curves) with 
Vp = 0.3. 



Figure 11: (Color online) 2Ao/A;bTc vs. AC/C„ for typical BCS superconductiors 
(squares) with the result calculated in the lA (the solid curve). 



Figure 12: (Color online) The temperature dependence of the magnetic moment ni 
calculated in the lA (solid curves) and FA (dashed curves) for q — 0.8 and q — 1.2: 
the result for g = 1.0 is plotted by the chain curve. 



Figure 13: (Color online) The q dependence of the Curie temperature Tq calculated 
by the lA (the solid curve) and FA (the dashed curve). 



Figure 14: (Color online) The temperature dependence of the electronic specific 
heat C calculated in the lA (solid curves) and FA (dashed curves) for q = 0.8 and 
q = 1.2: the result for g = 1.0 is plotted by the chain curve. 



Figure 15: (Color online) The temperature dependence of the inverse spin suscep- 
tibility 1/x calculated in the lA (solid curves) and FA (dashed curves) for q = 0.8 
and q — 1.2: the result for 5 = 1.0 is plotted by the chain curve. 



Figure 16: (Color online) The q dependence of the coefficients Cn,q for n = 1 — 4 in 
the generahzed Sommerfeld expansion calculated by the lA (sohd curves) and FA 
(dashed curves): note that ci,q = c^^q = in the lA whereas ci,q 7^ and cs^q 7^ in 
the FA. 



22 



This figure "figl.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig2.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig3.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig4.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig5.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig6.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig7.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig8.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "fig9.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "figl0.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "figll.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "figl2.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "figl3.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "figl4.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "figl5.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



This figure "figl6.gif" is available in "gif" format from: 



http://arxiv.org/ps/0906.0225v3 



